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Abstract 

Fully Developed Turbulence (FDT) occurs at the infinite extreme of the Reynolds 
spectrum. It is a theoretical phenomenon which can only be approximated experimen- 
tally or computationally, and thus its precise properties are only hypothetical, though 
widely accepted. It is considered to be a chaotic yet steady flow field, with self-similar 
fractalline features. A number of approximate models exist, often exploiting this self- 
similarity. We hereby present the exact mathematical model of Fractal Potential Flows, 
and link it philosophically to the phenomenon of FDT, building on its experimental 
characteristics. The model hinges on the recursive iteration of a fluid dynamical trans- 
fer operator. We show the existence of its unique attractor in an appropriate function 
space - called the invariant flow - which will serve as our model for the FDT flow field. 
Its sink singularities are shown to form an IFS fractal, resolving Mandelbrot's Conjec- 
ture. Meanwhile we present an isometric isomorphism between flows and probability 
measures, hinting at a wealth of future research.* 

MSC class: 76F20 (primary); 76F25, 37A30 (secondary) 

Keywords: turbulence, chaotic advection, attractors, fractals, invariant measure 
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1 Introduction 



We begin our discussion in Section 2.1 with an extensive overview of Fully Developed Tur- 
bulence (FDT), in an attempt to identify the definitive characteristics of this hypothetical 
phenomenon at the infinite end of the Reynolds spectrum, which can only be approximated 
and observed via experiments or simulation. We then proceed to give an account of the 
beautiful idea of Potential Flows and their superposition, heretofore seemingly unrelated to 
FDT, yet intriguingly sharing its inviscid and steady characteristics. We proceed to detailing 
the influential direction of non-deterministic experimental Fluid Dynamics, called Chaotic 
Advection, which we shall reformulate in the language of Optimal Transport Theory. Since 
the well-known IFS invariant measure is the fixed point of the iterative evolution under the 
measure transfer operator, it will serve as our inspiration for showing the unique existence 
of an invariant flow. Since the transfer operator distributes singularities over an IFS fractal, 
this adds impetus to pursuing a transfer model, considering Mandelbrot's Conjecture. In 
our remarks on this overview, we make a decisive step to discard Navier-Stokes evolution 
on mathematical grounds, observing the non-smooth nature of intermittency, and turn our 
attention to discrete inertial transfer evolution on piecewise steady flows. 

The fusion of the above theories will come naturally and effortlessly in the form of Fractal 
Potential Flows, introduced in Section 3. The random intermittent interaction of a finite 
system of eddies, is directly modelled via the recursive iteration of a weighted transfer op- 
erator in Section 3.2. The contractive action of each eddy sink on the flow field, is modelled 
with a corresponding pushforward transfer map, induced by a contraction map of an IFS. 
Since the collective action of the eddies is probabilistically weighted, it can be interpreted as 
an "expected value" expression, while physically speaking as the weighted effect each eddy 
has on the flow field. The iteration of the transfer operator, is reasoned to correspond to the 
non-smooth discrete energy bursts of intermittency. 

The evolution equation towards FDT having been reasoned to be the discrete transfer recur- 
sion, our goal consequently becomes to show that it converges to a unique attractor flow - the 
FDT flow field itself - in some appropriate function space. In finding the correct space, eddies 
remain our guiding inspiration, considering that a system of eddies becomes a superposed 
eddy when zoomed out, possessing a complex character at complex infinity as discussed in 
Section 3.3. In their local universe, the eddies interact to produce an increasingly fractalline 
flow field, approaching the desired invariant flow shown in Section 3.5. The required itera- 
tion number to reach FDT is infinite, while the total time can be potentially finite, if the 
time-spacing is for instance in geometric regression, as discussed in Section 3.6. 

Lastly, the sinks of the invariant flow field are shown to form an IFS fractal in Section 3.7.1, 
implying the ultimate conclusion that the geometrical study of such sets, is of fundamental 
relevance to analyzing an FDT flow field. 
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2 Preliminary Concepts 



2.1 Fully Developed Turbulence 

Turbulence is a challenging phenomenon to model, primarily because it is difficult to de- 
fine rigorously. One must isolate from experimental observation its definitive characteristics, 
before attempting to model it. However, even its main characteristics are difficult to iden- 
tify, and one may never feel secure that some crucial characteristic may have been excluded 
from the assumptions of a model, or that some assumption has been made which is not in 
agreement with reality. Furthermore, one must carefully differentiate between observed and 
intuitive characteristics as described by other researchers, as we hereby will. One must also 
restrict the problem of modelling sufficiently, in order for a specific solution to exist, which 
is why we shall focus on Fully Developed Turbulence (FDT). Although other approximate 
models may exist, one only needs to consider the phenomenon itself as defined from experi- 
ments, when attempting to build a new model. Currently FDT has no known exact model, 
which is a major issue considering its widespread occurrence in Nature and technology. 

A rich literature exists on the topic of FDT, and we recommend the overview of Ecke [6] 
for an introduction. Our overview assumes that the reader is familiar with elementary Fluid 
Dynamics. 

2.1.1 Definition 

Turbulence has been the interest of classical physics for centuries. Its name comes from the 
Italian word "turbolenza" given by Leonardo da Vinci, one of the first to study and depict 
the phenomenon. In terms of the Reynolds number, flows above an Re of 5000 or higher are 
considered to be turbulent. 

The phenomenon of FDT is an extremal case of turbulence, based on experimental obser- 
vations at high Reynolds numbers, typically above 10, 000. The words "fully developed" in 
FDT refer to its main theoretical characteristics: an infinite Reynolds number and a steady 
unchanging yet chaotic velocity field, independent of viscosity. This final unchanging nature 
of the field is meant to characterize an equilibrium state or invariance with respect to its evo- 
lution equation, whatever it may be. Furthermore, the velocity field, its set of singularities, 
as well as its evolution is hypothesized to be self-similar or fractalline in some sense. 

The precise ideal conditions for FDT to occur are currently unknown, and its exact causes, 
evolution, and characteristics are hypothetical. With some contemplation, one may realize 
that a mathematical theory is the only possible resolution. It is currently believed that the 
Navier-Stokes Equations can potentially describe turbulence at high Reynolds numbers. 
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2.1.2 Relevance 



The modelling of turbulence - whether approximate or exact - is of tremendous importance, 
considering its various and widespread applications. 

The study of atmospheric and oceanic turbulence occurs in the context of Geophysical Fluid 
Dynamics and Climate Theory - primarily in the form of 2D turbulence. Quantum turbulence 
in quantum fluids, such as superfluids, occurs at close to zero viscosity. Magnetohydrody- 
namic (MHD) turbulence is relevant to astro- and geophysics, ranging from the study of 
planetary magnetic fields to large-scale industrial facilities, and plasma flow in stars and fu- 
sion reactors. In biology, plankton distribution is influenced by turbulence, and is important 
since plankton comprise the bottom of the oceanic food chain. 

2.1.3 Characteristics, Causes, and Evolution 

Turbulence literature describes a set of observed and intuitive characteristics, and one must 
carefully distinguish between them. The causes and evolution towards Fully Developed 
Turbulence are also hypothetical, yet reasonable and experimentally inspired. Observed 
characteristics include high and erratic velocity fluctuation over a large range of coupled 
spatial and temporal scales. Consequently, a turbulent flow field exhibits structure at many 
length scales, often described as self-similar or scale-independent in some sense. 

The onset and evolution of turbulence can be fuelled by various circumstances, such as the 
interaction of eddies in the flow, by injected kinetic energy into the flow such as via stirring, or 
even by the shape or surface of the container. Examples include flow from a faucet, flow past 
rocks or a curved wall, or the rotation of the Earth. A periodicity in the circumstances often 
plays a role in inducing intermittency. Mathematically speaking, development can occur 
due to periodic changes in the external forcing, or variations in the boundary conditions. 
In intermittency, transitions are observed to occur abruptly to successively more complex 
states. The evolution of turbulence is suspected to be purely inertial, meaning due to the 
actual motion of the fluid, even if fuelled externally. Flows where inertial effects are small, 
tend to remain laminar. 

The induced eddies and their intermittent inertial interaction are considered to be the in- 
trinsic elements of turbulent dynamics. Since the injected energy at the top scale is periodic, 
the eddies in their evolutionary hierarchy receive periodic bursts of energy as well, causing 
their interaction to be intermittent. In between each energy burst, a brief quiescent period 
occurs, where the fluid flows along its natural course i.e. its steady streamlines. 

A hypothetical picture of the eddy hierarchy is the idea of vortex stretching and the Richard- 
son energy cascade. The eddies on the free surface of the fluid, will have axes roughly per- 
pendicular to the surface. A thinning of the eddies will occur due to the conservation of 
fluid elements. This results in the breakdown of larger structures into hierarchically smaller 
ones, causing a cascade of energy. This is a purely inertial and inviscid process, which con- 
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tinues until the local structures are small enough that viscosity causes their kinetic energy 
to dissipate as heat. 

Conjecture 2.1 (Mandelbrot [11]) At infinite Reynolds numbers (FDT), the energy dissi- 
pation of the fluid concentrates in a set of non-integer Hausdorff dimension. 

The conjecture has been resolved by V. Scheffer in 1980, though the proof is non-constructive. 
If one wishes to construct an FDT flow field, this conjecture may serve as an inspiration. 

The Richardson cascade is primarily a 2D phenomenological hypothesis, taking place on 
the surface of the fluid. Turbulent fluids typically occupy 3D space, but motion in the 
vertical direction can become less relevant, and we are left with a quasi-2D hydrodynamics. 
The degree to which the third dimension is suppressed, primarily depends on the ratio of 
the vertical to the lateral scales, and the effect of external forces causing dominant lateral 
velocities over vertical ones. 

2.1.4 Existing Models 

Previous efforts at modelling turbulence, and FDT in particular, have been approximate or 
statistical in nature, mostly in the form of scaling laws derived from the equations of motion. 
This direction was inspired by Kolmogorov's influential model based on the Richardson 
cascade, culminating in the current Beta Model developed by Novikov, Stewart, Mandelbrot, 
Kraichan, Frisch, Sulem and Nelkin et al. For a detailed overview, we recommend Jou [8]. 
There is also the computational direction of simulating large Reynolds number flows on 
high-performance computers, originally suggested by John von Neumann. 

The theory of Chaotic Advection is an experimental and mathematical direction, meant to 
investigate the intermittent interaction of eddies directly. Its language is Potential Flow 
theory and Dynamical Systems. We will dedicate an entire section to its review. 

Detailed understanding and prediction from first principles still elude turbulence theory, 
which mostly bases itself on the Navier-Stokes Equations. Currently, it is considered to 
be highly unlikely that an exact rigorous mathematical model should be possible for either 
turbulence and its evolution, or its end result, Fully Developed Turbulence itself. If such a 
model exists, it is conjectured to be geometrical. Nevertheless, an exact model must only be 
in agreement with the observed phenomenon itself, and not with other derived models. 

2.2 Potential Flows 

For laminar flows outside the boundary layer, potential flows are considered to provide a 
sufficient model. Yet the turbulent boundary layer itself is often modelled approximately 
with a superposition of vortices (Vortex Dynamics). These flows are assumed to be steady 
(time- independent), ideal (zero viscosity / inviscid, uniform density, and incompressible), 
and irrotational (Lebesgue almost everywhere, denoted as a.e.). 
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By these conditions, conservation of mass for a velocity field v = (v 1: v 2 ) : M 2 -> M 2 is 



= div i> = + (92t>2 a.e. 



and the equations of motion become 

1 , ,o . 1 

-V u — v x curl v = — vi 

2 11 Q 

where g is the density and p is pressure, further reducing to 




under the irrotationality requirement 

= curl v = d 2 v 1 — d\V<i a.e. 



Supposing that ip, <p G C 2 e (M 2 , M) are harmonic conjugates, meaning they satisfy the Cauchy- 
Riemann Equations 

Aip = A(f) = a.e. and di4> = d 2 vp, d 2 <p = —diip a.e. 

then v := — V0 = lVip satisfies mass conservation and irrotationality a.e. The pressure field 
can be calculated from v as above. We call such a.e. harmonic ip f° r which a conjugate 
exists, the stream function of a potential flow, and <fi the potential function. We note that 
the potential function corresponding to a stream function is only unique up to a gradient. 
Nevertheless we denote it as :— (f), and denote equivalence in the gradient as 

0! = 02 V0i = V0 2 a.e. 

A harmonic function is known to admit a conjugate if its domain is simply connected. 
Furthermore ifj = —ip. 

By their relation to the velocity field, the curves of constant vp represent the streamlines of 
the flow, and the curves of constant the equipotential lines. Having two pairs of harmonic 
conjugates ipi,(fii and ip 2 ,(p 2 their linear combinations aipi + bip 2 and a<p\ + b(fi 2 for any 
a, b G K. are also harmonic conjugates. So the principle of superposition holds, as long as the 
boundary conditions are also correspondingly combined. Since the equations of motion have 
a unique solution for a set of boundary conditions, and they translate to Laplace's equation 
(the solution of which for appropriate boundary conditions also exists and is unique up to a 
gradient), we may conclude that potential flows fully characterize the set of all flows which 
are steady, ideal, and irrotational a.e. 

Elementary potential flows, which are often superposed to create more complex ones, include 
sinks, sources and the (circular) vortex. The stream and potential functions of a source 
(q > 0) or a sink (q < 0) with strength q G K. are defined at z G C as 

q q 
ip(z) = — — Arg z mod \q\ and <p(z) = —— In \z\ 

2,71 ' 2,71 
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while for the orthogonal circular vortex 

r r 

ip(z) = — In |2;| and <p(z) = — — Arg z mod |T| 

where T is the circulation around any closed path containing the vortex (strength of the 
vortex) . 

Superposing a sink and a circular vortex with a general center p G C results in a logarithmic 
vortex (eddy) as follows 

ip(z) = — — Arg(z — p) + — In \z — p\ mod \q\ 

2,71 2,71 

(f)(z) = — ^- In \z — p\ — —Axg{z — p) mod |T| 
2tx 2ix 

Since ip is differentiable in the M 2 sense almost everywhere (except at p and along the branch 
■ip(z) = q) we have that the corresponding velocity field (extended continuously to the entire 
plane) is the following 

Z7T \Z — p\ \Z — p\ 

Here we consider the gradient vectors to be on the complex plane and differentiation in 
the bivariate sense, meaning = diip + id 2 ip. In general, we will sloppily identify z 
ip(z), z G C with (x, y) ip(x + yi), (x, y) G M. 2 when it is more convenient to do so. 



2.3 Chaotic Advection 

2.3.1 Aref's Blinking Vortex-Sink System 

Advection is the idea of a fluid transporting light matter (typically tracer particles) on its 
free surface, so that the velocity of the particles is given by the velocity field of the flow. 
If if) is the stream function of the velocity field, then this corresponds to the Hamiltonian 
equations v\ = — c^, v 2 = ditp. Hence the configuration space of an advected particle is 
the phase space of this Hamiltonian system. Often such a system exhibits chaos, or sensitive 
dependence on initial conditions. The notion that laminar flows can produce chaotic particle 
trajectories, is considered today to be a cornerstone of Fluid Dynamics. Hassan Aref is the 
originator of this idea, and his influence has propagated far. 

The flows mostly studied in this context, are piecewise steady flows. Meaning that on the 
evolution timeline of the flow, the velocity field is steady for disjoint consecutive time in- 
tervals, and a non-smooth abrupt transition (jump) to a new velocity field occurs at each 
partitioning moment in time. The jumps act along the streamlines of two circular (or log- 
arithmic) vortices, which take alternating (non-random) turns in perturbing the flow field. 
This is called the Blinking Vortex(-Sink) System, and has been examined extensively by a 
number of theorists and experimentalists. It is meant to directly simulate the intermittent 



8 



interaction of (circular or) logarithmic vortices, which is considered to be the elementary 
driving mechanism of the evolution towards Fully Developed Turbulence. 

Aref in his highly influential paper [1] originally defined the setup with two circular vortices, 
which has been modified in [2] to two vortex-sinks (eddies), and was further examined by 
Karolyi and Tel [10] as well as a number of other enlightened researchers. 

In particular, Wiggins [20] discusses the method and implications of modelling this blinking 
as the alternating (non-random) infinite iterative composition of two maps (called linked 
twist maps, the study of which goes back to Devaney [5]) in order to show the existence of a 
chaotic invariant set via the Conley-Moser conditions. These ideas are further examined with 
Ottino in [21], and bear an inspirational resemblance to our own, to be presented here. For 
an overview of Chaotic Advection (or Mixing), we recommend the book by Ottino [13]. 

2.3.2 Ott's Chaos Game 

Ott et al. [22] has made a significant improvement on Aref 's model. The regular iterative 
alternation between two sinks has been modified to a random iteration over a finite family of 
mappings that transform the flow at the intermittent time partitions, in between steady flow 
periods. So at step I G N (time partition ti) the state of the flow o\ (usually the position of a 
tracer particle on the fluid surface) is transformed as 07+1 = Tj(o'j) where T| G {Ti, . . . ,T n } 
is chosen at random, and where the precise form of T; can vary and is often left implicit. 
This random discrete evolution is reminiscent of the Chaos Game of Barnsley [3] . Though as 
we will see, this formulation can be closely aligned with experiments, it is difficult however 
to treat it mathematically. Ott et al. further assume that t\ are spaced evenly and their 
stream functions are time-dependent, as typical in the theory of Chaotic Advection. However 
as we will see, both of these are unnecessary constraints, in our explicit and more natural 
formulation of this chaos game, in Section 3.2. An overview of the work of Ott and Sommerer 
et al. in collaboration with Tel has been presented in [16], building upon [22, 15, 14]. 

2.3.3 Sommerer's Experiment 

Chaotic Advection is a de facto experimental theory, and its predictions have been verified 
by both computational and real-world experiments. Sommerer's apparatus [15, 16] carried 
out at the JHU Applied Physics Laboratory, consists of a tank of water with a set of sources 
at the bottom, which transmit pulsing injections from an external pulsatile flow modulator 
that randomly regulates which source receives the next injection at equal time intervals, 
circulated by an external pump. Each pulse corresponds to an iterative step 07+1 = TJ(cr/) 
where the fixed set of TJ G {Ti, . . . , T n } represents the actions of each source on the particles. 
These fluorescent particles on the fluid surface are then traced over their chaotic paths. In 
a sense, this apparatus is an inversion of the Aref setup, since it utilizes sources instead of 
sinks. Nevertheless, its relevance lies in attempting to simulate the intermittent evolution of 
turbulence, in the spirit of Ott's chaos game recursion. 
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2.4 IFS Fractals 



The attractors of Iterated Function Systems - IFS Fractals - were pioneered by Hutchinson 
[7], further discussed by Barnsley and Demko [4], and may be the most elementary fractals 
possible. They are the attractors of a finite set of affine linear contraction mappings on the 
plane - the "function system" - which when combined and iterated to infinity, converges to 
an attracting limit set, the IFS fractal itself. Our formulation follows [17]. 

2.4.1 Definition and Existence 

Definition 2.1 Let a 2- dimensional affine contractive mapping (briefly contraction or con- 
traction map) T : C -> C be defined for all z G C as T(z) := p + ip(z — p) where p G C is the 
fixed point of T , and <p = Xe^ 1 G C is the factor ofT, with X G (0, 1) the contraction factor 
of T , and $ G (— tt,tt] the rotation angle ofT. 

Note that an equivalent definition may be given using unitary rotation matrices R G M 2x2 , 
R T R = I corresponding to e . Then the contraction maps take the form 

T(z) =p + XR(z -p) (z G R 2 , p G IR 2 , AG (0, 1)) 

This version shall be useful for certain proofs in later sections. 

Definition 2.2 Let a 2-dimensional affine contractive n-map iterated function system (briefly 
IFS or n-map IFS, n G N) be defined as a finite set of contractions, and denoted as 
T := {Ti, . . . ,T n }. We will denote the set of indices as Af := {1, . . . ,n}, the set of fixed 
points as V := {p±, . . . ,p n }, and the set of factors as $ := {tpi, . . . , <£>„}. 

Definition 2.3 Let T = {Ti, . . . ,T n },n G N be an IFS. We define the Hutchinson operator 
H belonging to T as 

n 

H{S) := |J T k (S) where T k (S) := {T k {z) : z G S} for any S C C 

k=l 

and call H(S) the Hutchinson of the set S. 

Theorem 2.1 (Hutchinson [7]) For any IFS with Hutchinson operator H, there exists a 
unique compact set F C C such that H(F) = F. Furthermore, for any compact S C C, the 
recursive iteration S n+ \ := H(S n ) converges to F in the Hausdorff metric. 

Proof The proof follows from the Banach Fixed Point Theorem, once we show that H is 
contractive in the Hausdorff metric over compact sets, which form a complete space. □ 

Definition 2.4 Let the set F in the above theorem be called a fractal generated by an IFS 
with Hutchinson operator H (briefly IFS fractal). Denote (Ti, . . . ,T n ) = (T) := F. 
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2.4.2 The Address Set 



The address set arises from the iteration of the Hutchinson operator, and it is a way to label 
each fractal point. Since we can start the iteration towards F with any compact set, we will 
choose the primary fixed point for simplicity, which is any point in V of our preference. 

Definition 2.5 Let M L :— M x . . . x Af be the index set to the L-th Cartesian power, and 
call this L E N the iteration level. Then define the address set as A := {0} U UzLi N L . For 
any a E A denote its k-th coordinate as a(k), k E N. Let its dimension or length be denoted 
as \a\ E N so that a E A/"' a ' and let |0| := 0. Define the map with address a E A acting on 
any z E C as the function composition T a (z) := T a m o . . . o T a n a ^(z). Let the identity map 
be T := Id. For a set of weights {w\, . . . , w n } C (0, 1) define w a := w a (i) ■ . . . ■ w a (\ a \), for a 
set of factors {<p u . . .,tp n } define tp a := v? a (i) • • • ■■ ( fa{\a\), and for a set of angles {A, . . . ,# n } 
define $ a := $ a{l) + ... + i? a (|a|)- 

Theorem 2.2 For any primary fixed point p EV we have 

F = (T 1; . . . , T n ) = lim H L ({p}) = C\{T a (p) : a E A} = C\{T a (p k ) : a E A, p k EV} 

L— >oo 

We call this the address generation of F. 

Proof The proof follows from Theorem 2.1 with either of the initial sets {p} or V. □ 

Definition 2.6 Inspired by the above theorem, we will call F L := {T a {p) : a E A, \a\ = L} 
for some L E N and primary fixed point p EV the L-th level iterate towards (T 1; . . . , T n ) : = 
{T a (p) : a E A} the countably infinite IF S fractal. 
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2.5 Optimal Transport over Measures 



Transportation Theory, often referred to as the theory of Optimal Transport, is the study 
of the optimal transportation and allocation of resources. Given two finite sets of producers 
and consumers equal in number, with each producer supplying a specific consumer according 
to a bivariate cost function for transportation, the goal is to find the minimal arrangement 
between producers and consumers, called the optimal transport plan or map. 

The discrete problem can be generalized to probability measures - the Monge-Kantorovich 
formulation, named after its two most influential creators. Like many optimization problems, 
this minimization problem has a dual maximization equivalent. When the cost function 
over a compact metric space (X, dx) is the metric itself, then the dual becomes the first 
Wasserstein distance (also called Hutchinson metric) between the elements of the set of 
Borel probability measures G M-x over X, defined as 

d{v x , v 2 ) = sup I jf / dfa -u 2 ):fe C X {X) | 

where C\{X) are the Lipschitz-1 functions over X. That this is indeed the dual to the 
minimization problem, was shown by Kantorovich and Rubinstein [9]. Furthermore (Aix, d) 
is known to be a complete metric space. 

Generally the transport plan is a pushforward map (transfer operator) of the form TV = 
v o T _1 [y G M-x)- If T* is contractive in d, then by the Banach Fixed Point Theorem it has 
a unique invariant measure V = T*T>. Taking a finite weighted average of such pushforward 
operators, the resulting transfer operator 

T = u>iT* + . . . + w n T* where w k G (0, 1), W\ + . . . + w n — 1 

will also be contractive in d over Aix, and possesses a fixed point. Thus the following 
theorems hold for an IFS {Ti, . . . ,T n } in the form T k (z) = Pk + ^kRk(z — Pk), z E M. 2 
supposing that F C X C 1R 2 . 

Theorem 2.3 (Hutchinson [7]) For any IFS {Ti, . . . ,T n } and weights {wi, . . . ,w n }, the 
transfer operator T = W\T^ + . . . + w n T* is contractive on (Aix,d) and attains a unique 
invariant probability measure V = Tu over the compact set X C 1R 2 ; with support (Ti, . . . , T n ) . 
We call this the invariant measure with respect to the IFS {T 1; . . . ,T n }. For any u G M.x 
the recursion vl+i = wiT^vl + . . . + w n T^v^ G Aix, L G N converges to the invariant 
measure. 

We remark that with the primary fixed point p and the notations of Section 2.4.2, the 
measures 

= S p , vl=Y^ w " T a 5 P e M x, LeN 

\a\=L 

satisfy the above recursion, and so v{S) = lim(i/^) as discussed in [17]. We also observe that 
the support of v>l is concentrated on the L-th iterate Fl := H L ({p}) = {T a (p) : a G A, \a\ — 
L} approaching the compact IFS fractal F = (Ti, . . . , T n ). 
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We further remark that defining 

V x := {g G C ae (R 2 ,R) :g>0, Supp g C X, [ g = 1} 

Jx 

V x ■= Y e M x ■ 3g G V x V5 C X Borel set : u(S) = J g J 

and keeping the above metric d, then Conv{5 p : p G IR 2 } C Vx C Aix and V G Px, so it 
remains true that there exists a unique invariant measure in (V x ,d), and it is the same as 
the earlier. Furthermore T : Vx Vx since Supp(TV) = //(Supp u) C X if v G "Px- 

The study of transfer operators - also called Ruelle (-Perron- Frobenius-Markov) operators - 
is a rich field, and their largest eigenvalue is typically one, while their eigenfunctions are 
usually fractalline or self-similar in some sense. This has profound implications for classical 
mechanics, such as the increase of entropy or the irreversibility of time. An eigenvalue of 
one, corresponds to a state of equilibrium. 

For a more complete discussion of Optimal Transport, we refer the reader to Villani [18] and 
McCann [12]. For an overview of IFS, we recommend articles [19, 7, 4]. 

2.6 Remarks 

We hereby wish to make some remarks on the previous sections, which will shape our ap- 
proach to be presented in Section 3. 

First of all, it seems reasonable to raise the doubt whether Fully Developed Turbulence 
in full three dimensions even exists. Certainly, flows can approach above 10, 000 Reynolds 
numbers, but it is doubtable whether the other defining characteristics of FDT, including 
the equilibrium state of an unchanging velocity field at Re = oo that is independent of 
viscosity, can evolve in a fully three-dimensional flow field. Furthermore, we recall that the 
Richardson cascade that inspired the notion of self-similarity, is primarily a 2D picture on 
the free surface of the fluid. It seems that 3D space has one too many degrees of freedom, 
while a limiting 2D free surface forces the fluid to self-organize via self-similarity. Thus we 
restrict our attention to the 2D free surface in our effort to model FDT, in the spirit of 
Chaotic Advection. 

Secondly, we point out the mathematically unreasonable effort of various other approaches, 
which attempt to model FDT via the Navier-Stokes Equations. The discussed energy bursts 
fuelling the erratic characteristic of intermittency, are clear signs of the non-smooth na- 
ture of the evolution towards FDT. The Navier-Stokes Equations however, require temporal 
differentiability of the velocity field, which does not hold by the observed phenomenon of 
intermittent evolution, summarized above. Therefore one has no choice mathematically, but 
to discard the Navier-Stokes Equations, in hope of a more fitting evolution equation, in line 
with the experimental picture of intermittent inertial evolution. 
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Thirdly, we point out that prior to the advent of Chaotic Advection, potential flows had 
been considered to be useless for the study of turbulence. Yet ironically, potential flows 
share with FDT the characteristics of independence with respect to viscosity and time, 
both being steady flow fields. In fact for increasing Reynolds numbers, the Navier-Stokes 
Equations approach the inviscid Euler Equations, implying that the flow would converge to 
the inviscid solutions of Potential Flow Theory. Thus the Reynolds spectrum appears to be 
an ouroboros, similarly to the Riemann circle. 

Lastly, we point out that in the discussion of Chaotic Advection, we have not introduced the 
standard mathematics in full, generally employed in the theory. This is partly because we find 
the existing methods to be suboptimal for actualizing the potential in its concepts, and partly 
because the theory in its current form is often intuitive and in need of a rigorous mathematical 
language. So we feel that the theory needs a complete reformulation unifying the Aref and 
Ott directions in the spirit of Wiggins, which is one of our goals in this paper. 

So to project our intentions, the model to be presented shall employ the methods of Optimal 
Transport Theory to reformulate in a new language the theory of Chaotic Advection, in the 
setting of Potential Flows, with the purpose of creating an exact geometrical model for the 
evolution towards the final steady equilibrium state of Fully Developed Turbulence. 
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3 Fractal Potential Flows 



We have introduced above the theoretical elements that we intend to weave together in this 
section and the next, in our effort to create an exact model for Fully Developed Turbulence 
and the evolution towards it. We will introduce Fractal Potential Flows as our exact model, 
which arise due to the intermittent transfer interaction of eddies. We define a weighted 
transfer operator over a certain class of stream functions, and show the existence of its 
unique fixed point - the invariant flow - characterized by a stream function whose rotated 
gradient field will serve as our model. The set of sink singularities of this flow field is shown 
to be an IFS fractal, explicitly resolving Mandelbrot's Conjecture. 



3.1 Eddy Invariance 

In Section 2.2 we have defined the stream function of an eddy with sink strength q = Chi A 
and vortex strength or circulation F = Cd as 

ip(z) = Arg(z — p) -\ m|2 — p\ mod —Gin A 

2n 2% 

where <p G C, A = \<f\ < 1, $ = Arg (p and C > is an arbitrary parameter. The associated 
velocity field (continuously extended to the entire plane) becomes 

GLog <f 1 z — p 



2n \z — p\ \z — p\ 



We observe that the eddy stream function is invariant under the contraction T(z) = p + 
(p(z — p), meaning T*ip = if). This is expected geometrically, since the orbit T 1 (zq) = 
p + ^(zq — p), t G R traces out a logarithmic spiral centered at p G C, and the streamlines 
of constant if) are spirals of the same pitch. So in other words, if> is the fixed point of the 
pushforward transfer operator T*. This raises the question whether if> is unique, and if not, 
then what space would guarantee its uniqueness. Defining ipo(z) := Arg [z — p) mod $ we 
see that T*ip Q = ip also holds, therefore ip cannot be unique in such a general setting. 
So our ultimate goal becomes finding the proper function space where uniqueness can be 
guaranteed, with a reasonable physical interpretation. We begin our search by deriving 
some fundamental properties of transfer, which will have profound implications. 

Theorem 3.1 For a similarity contraction of the form T(z) = p + (p(z — p), \f \ < 1 the 
following properties hold (with differentiation in the 1R 2 sense, and v[ip] being the continuous 
extension ofiVip). 



AT* = ^ T*A and T*ip = T*4> and v[T*if)} = ^ T*v[ 



Proof We show the above in the R sense first, which translates to the complex sense, 
using T(z) = p + XR(z — p), p,z6l 2 where R is the rotation matrix corresponding to e" 01 . 
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We first show the third property, keeping in mind that v [^K- 2 ) — iVip(z) a.e. Gradient is a 
column vector in IR 2 , which corresponds to the complex vector d\ + 182- Differentiation D 
results in a row vector however, so we must take a transpose, meaning V = D T . Applying 
the generalized chain rule, we get 

VO o T- 1 ) = DO o T~ X ) T = ^(D^ o T- 1 ) ■ \ R ~ l ^J =^R-Vtpo T _1 

Multiplying by a 90° rotation (or i), we get the third property. To get the first property, we 
observe that A = tr DV, so by the chain rule, the properties of trace, and the above 

AT*ij = tr DVT*V = tr D f^R ■ T*V^J = tr Qi2 • T*DV^ • = ^ T * A ^ 

The second property also follows from the third, observing that in the above complex sense 
ip and if) are conjugates iff = iV^ by Section 2.2. So taking T* of both sides, multiplying 
by i and dividing by A, we get that V(T*ip) = iV(T*^) implying that V(T*^) = V(T*ip) 
and thus T*ip = T*^p. □ 

Among many things, the first property also implies that if ip is a.e. harmonic then so is T*ip. 
Together with the second property, this implies that the space of potential flows is closed 
under pushforward transfer by a similarity contraction. 

We end our discussion of a single eddy by emphasizing the underlying idea that one applica- 
tion of T* to a flow field represented by the stream function, corresponds to one intermittent 
energy burst fuelling an eddy to transform the flow field along the logarithmic spiral orbits 
of T. If T* can be shown to be contractive over stream functions in some complete space, 
then the final equilibrium fixed point stream can also be shown. 



3.2 The Transfer Operator 

Considering that potential flows form a vector space as discussed in Section 2.2 and by 
Theorem 3.1, the weighted transfer operator 

Tip = WtT^tp + . . . + w n T*ip 

for any IFS {Ti,...,T n }, T k (z) = p k + (p k {z - p k ) and weights w k E (0, l),^ k w k = 1, 
maps the stream function of a potential flow to another such flow. One application of T to a 
flow field corresponds to the random application of some Tj* e {Tj*, . . . , T*} with probability 
P(Ttp(zo) = T^ip(z )) = w k at some point z in the plane. Indeed, the iteration of the 
transfer operator T corresponds to the intermittent randomly alternating contractive action 
of a system of n eddies along the orbits of T k . Here the flow field is piecewise steady in 
between the energy bursts fuelling some particular eddy at a certain intermittent moment. 
During the steady periods, the flow field is of the form T L ipo, L G N for some initial flow 
■00- The partitioning intermittent moments in time ti can be spaced in any way over the 
timeline. This constitutes our mathematical reformulation of Chaotic Advection. 
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We see that after the L-th iteration we have 

\a\=L 

where T a and w a were defined in Definition 2.5. Probabilistically speaking, this implies that 
ipij is the stochastic superposition of n L eddies with transfers T*, transforming the flow field 
at the L-th intermittent moment with probabilities P(iPl(z ) = T*if) (z )) = w a . We also 
note that by Theorem 3.1 the corresponding velocity transfer operator is 

v [T^j] = w 1 y^ T?v[ip} + ... + w n r^^ T*v[ip] 

There are two possible interpretations of the iteration of the above transfer operator - prob- 
abilistic and physical. Probabilistically it means that considering the next iterate ipL+i to be 
a random variable, its expected value is TipL. A related financial interpretation is that the 
collective action of the eddies, is a portfolio. Physically, weighted transfer can be interpreted 
as weighting the effect each eddy can have on the fluid at an intermittent moment. 

Succinctly, we observe that iteration corresponds to intermittency, and hereby theorize that 
the intermittent evolution of turbulence towards the equilibrium state of Fully Developed 
Turbulence, corresponds to the iteration of T ad infinitum to an invariant flow Tip = ip. 
Therefore ip will represent our model for the fully developed turbulent flow field, and show- 
ing its unique existence now becomes our primary goal, along with finding an appropriate 
function space where this is possible and physically most reasonable. 

Theorem 3.2 For the weighted transfer operator T the following properties hold. 

n n 

AT = Yl 5 T k A and ^ = T ^ and v l T ^ = Yl Wk r% T * v ^ 
k=\ k k=i 

Proof The proof follows trivially from Theorem 3.1. Note that the first property implies 
the preservation of harmonicity a.e. or physicality upon transfer, as discussed in Section 2.2. 
□ 

3.3 Flow Character 

In constructing the proper flow space, our inspiring objective is to ensure the unique existence 
of an eddy as the attractor of its generating pushforward transfer map. In resolving this 
question, we show a correspondence between eddies and the Dirac delta function, which hints 
at a possible general correspondence between stream functions and the density functions of 
probability measures, via Poisson's Equation. Finding the ideal kind of boundary condition 
will prove crucial to our quest, and meanwhile the proper space of flows shall slowly reveal 
itself to us. 
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First of all, we observe that the velocity field of the eddy discussed in Section 3.1 is char- 
acterized by the complex parameter c = |^Log cp which can be extracted from the stream 
function ip by the operation 

v[ib](z) z — p = c -, ~io z — p = c 

\z — p\ l 

We will denote the stream function of an eddy with character c G C, Re c < and centered 
at p G C, with q c := — 27rRe c, T c := 27rlm c as 

q c p c 
AA Z ) '■= w~^i. z - P) + m \ z ~ P\ mod Qc and ip c := ip cfi 

Certain flows may have a similar "eddy character" when their flow field is zoomed out, even if 
they exhibit varying streamline behaviour locally around the origin, as defined below. 

Definition 3.1 We say that a stream function ip : C ->■ M satisfies the boundary condition 
at infinity with character c G C, if for any p G C and any sequence (zj) C C, \zj — p\ -> oo we 
/iai>e 3 lim^oo t> [^>] 2; j — p = c where the limit is taken in the complex sense. We denote 
this property as ip G BC™ or char(-?/>) = c. 

We remark that character is independent of the choice of p G C. To see this, take 



lim = lim 

j^OO j-iOO 



Zj Zj - p 



z,-p 



lim 1 = 

i^oo \Zj — p\ 



so lim(f [V'K^j)) = 0. Taking some other p' G C we have for the above (zj) that \zj —p'\ ^00 
since \zj — p\ < \zj — p'\ + |p' — p| and that 



Zj) Zj — p' — v[ijj](zj) Zj — p + v[ip](zj) p — p' -¥ c + = c as j ->• 00 

Theorem 3.3 Character is a linear map that is invariant under weighted transfer, or con- 
volution with a density function over the plane. 

Proof The linearity of char is trivial by definition. Let us suppose that char(?/>) = c. Since 
for some p and (zj), \zj — p\ -> 00 sequence Zj —p = ifk(T k ~ 1 (zj) —pk) + (pk —p) and applying 
Theorem 3.1 we have that 

n 

<fk 

k=l 



v[T^}(zj) Zj-p = v[ip](T k 1 (z j )) ip k (T k 1 (zj)-p k ) + (p k -p) 



J2 Wk J^ v ^ T k \ z i)) T k \ z j) -Pk + J2 Wk U^ v Wi( T k 1 ( z i)) Pk-P 



k=X ^ S Wk\ 



which approaches J2 k w^c + J2 k = c as j -> 00 meaning char(T0) = c. 

Now let q be a density function over M 2 and (zj) be the above sequence. We employ scalar 

and complex vector Riemann integrals as follows. 
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+ 



V('0 * q)(z) — / Vip(z — w)g(w) dw 
Jm 2 

v[ip * g](zj) Zj — p = / iVip(zj — w) (zj — w) — p g(w) dw + 
Jr 2 

/ iVvp(zj — w) w g(w) dw -> / cg(w) dw + 

JR2 JR 2 



c as j 



oo 



□ 



In what follows, we will denote with L°°(IR 2 ) the set of bounded functions over the plane, 
and with C 2 e (IR 2 ) the a.e. defined and twice a.e. continuously differentiable functions ip, for 
which iVip can be continuously extended to the entire plane, denoted as v[fa, furthermore 
the set of singularities of v[fa\ is compact and nonempty. 

Theorem 3.4 For any density function g over the plane and any c G C, Im c ^ 0, there 
exists aip G L°°(IR 2 ) flC 2 e (IR 2 ) nBC^° for which Aip = T c g a.e. Furthermore, this ip is unique 
up to a gradient (&.e.) and ip = fa * g. 

Proof The existence of such a ip is guaranteed using the Green function method of convolv- 
ing the fundamental solution fa (Aip c = T c 5q a.e.) with g which will also have a character 
of c by the previous theorem. So we turn to the question of uniqueness (up to a gradient). 
Let us suppose indirectly that 3^1,2 £ BC^°,ipi ^ fa such that Afa = Y c g = Afa a.e. Then 
Aip = a.e. with ip := fa — fa. By the Divergence Theorem, for some compact S C IR 2 with 
piecewise smooth boundary 

\Vfa 2 = [ |W>| 2 + M^ = / div(^W>) = / ^(V^,u) 

IS JS JS JdS 

Since fa >2 are bounded, then so is ip with say \fa < M. So for any sequence of disks 
Sj := B(p, rj), < r j < Tj + i, j G N, lim(rj) = 00 for which S^ C Si, dS^ D dSi = 



[ IWI 2 = 


[ ip(Vi>,u) 


<M 1 


/Sj 


JdSj 


JdSj 



Vfaz), 



p 



\z — p\ 



dz< M 



\v[fa(z)\ dz 



dSi 



where the integrals are in the usual IR 2 sense. If we can show that this last integral vanishes 
as j -> 00, then J R2 \Vip\ 2 = and so Vip = a.e. which would contradict fa ^ fa. 
Let us define Zj := axgm&x z€ds \v [fa(z)\, j G N which exist by the Extreme Value Theorem, 
considering the compactness of dSj and the continuity of v[fa on dSj C IR 2 \ S^. Since 
Zj G dSj, j G N we have \zj — p\ — Tj -> 00, j -> 00 which by char('0i i2 ) = c implies that 

1 

\v[fa( Zj )\ dz = \v[fa( Zj ) Zj -p\ 



v[fa(z 



3! z j 



p -> so 

v[fa(z)\ dz < 



OS, 



dSi 



\Zj-p\ 



2irrj -> as j -> 00 □ 



By the above results, we may conclude among many things that the eddy fa jP uniquely cor- 
responds via Poisson's Equation to the Dirac delta function 8 P among all bounded functions 
of character c, where the uniqueness is guaranteed up to a gradient. 
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3.4 The Flow Space and the Measure Map 

Definition 3.2 Let X C M 2 be a compact set such that F = (Ti, . . . ,T n ) C X and c G C 
some number for which Re c < 0, Im c ^ 0. For a function ip G C 2 e (IR 2 ) define 

M>](S) = fi^KS) := ±- I (5 CM 2 ) 

Let the set of flows with character c be defined as 

^>x,c := {ip e C(R 2 ) n : Ran V C [0, g c ), /#] G Px} 

VFe wi// ca// /x : ^ x,c Vx the measure map, and fjL[ip] the measure induced by ip. 
Recalling the metric d from Section 2. 5, for any fa^ G ^x,c denote 

D(fa,fa) = D X)C {fa,fa) := d{^[fa]^[fa]) = sup Ur Jj A(Vi - V> 2 ) : / e £i(X)| 

^4s m Section 2.2, we will consider two functions to be congruent fa = ip 2 iff^ipi = Vfa a.e. 
over X. Let the equivalence class of ip be denoted as [ip] := {fa £ ^x,c '■ 4> = ^0} an d let 

[^x,c] ■= M : i> e *x,c} = *x,c/Ker(fi c ) 

Further denote D([fa], [fa]) = ^([^i], [fa]) ■= D x , c {fa,fa). 
Then we call ([^x,c]i D) the flow space over X with character c. 

Theorem 3.5 ({$/ x ,c],D) is a well-defined metric space. The measure map fi : [^x,c] Vx 
is an isometric isomorphism, and the weighted transfer maps T : ^x,c ^ x,c- Furthermore, 
the two maps commute. 

Proof First of all, we observe that /1 : [^x,c] "Px is clearly an isometry, since by definition 
d(ji[fa],n[fa]) = D([fa],[fa]). H is bijective iff G V x 3![V] e : /#] = i>- The 

equality n[ip] = v here means that with the density function g belonging to v G Vx 

MKS) = p / = = /" e V5cl 

c J S J S 

which is equivalent to Aip = T c g a.e. on X. The unique existence of such a t/> G ^x,c f° r anv 
was shown in Theorem 3.4 up to a gradient, meaning up to a congruence or equivalence 
class. So we have that /x is a bijective isometry, also called an isometric isomorphism. 
For the well-definition of the space, we must show that D is a metric over the set 
D is positive, symmetric, and inherits the triangle inequality from d. The only question 
that remains is whether D([fa], [fa]) — implies [fa] = [fa]. Since = D([fa], [fa]) = 
d(fi[fa], fx[fa]) we have that /i[fa] = [J>[fa], and since \i has been shown to be bijective, this 
implies that [fa] = [fa] or fa = fa. 
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We go on to showing that T and fi commute, meaning T/x = fiT. We recall from Theorem 
3.2 that 

n 

A(T0)(z) = £^ {^)oT k \z) 

k=l A k 

and also see that for any S C X C M 2 and k G {1, . . . , n}, by taking Riemann integrals over 
M. 2 we have that 

nim(S) = 1 / A^oT, 1 ) = 1 / ^(AVO^C*)) dz = 

1 c l c Js A k 

= 1 / A^H d^ = /^IC^)) = (T fc >[^)(S) 

so taking a convex combination and using the linearity of \x we have the desired property 
/i[T-^] = T/i[^], ip G *&x,c or succinctly /iT = T/x. 

Lastly, we show that T maps from ^x,c to itself. T preserves C 2 e {R 2 ) since = H(S^), 
and it also preserves character according to Theorem 3.3. Since T is a convex combination, 
it preserves that Ran(T0) C [0, q c ). Furthermore, we have that //[T^] = T/i[-?/>] G Vx since 
G Vx and the measure transfer is known to map T : Vx -* Vx by Section 2.5. □ 

Physically speaking, the above space ([^x.cL D) represents the fact that the intermittent 
interaction of a system of n eddies, can be considered in a locally isolated manner over a 
compact X C M 2 of reasonable scale, [^x.c] versus *&x,c signifies that it is only the velocity 
field induced by a stream function that matters. Prior to the start of their interaction, the 
eddies ip Ck ,p k (notation of Section 3.3) superpose initially in a flow with a certain character c 
as follows. Letting ip := J2k w kTk?Pc k ,p k then by Theorem 3.1 

n n 

so clearly c = char(^> ) = J2k w k c k- As the intermittent interaction of the eddies progresses 
according to the transfer operator T, their overall character will be preserved to be this c by 
Theorem 3.3. 
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3.5 The Invariant Flow 



Theorem 3.6 The transfer operator T is contractive over ([^x,c], D). For any if>o E ^ x,c 
the iteration ipi+i = TipL G ^x,c converges to 

•0 = i\) c * V = lim } W a lf) c {- ~ T a {p)) E ^x,c 
\a\=L 

for which Tip = if), and here V is the invariant measure of Theorem 2.3 and p is the primary 
fixed point. We will call this if) the invariant flow (or fractal potential flow) with respect to 
the IFS {Ti, . . . , T n } and character c E C. 

Proof The contractivity of transfer over stream functions, follows directly from that of the 
measure transfer in Theorem 2.3, considering that the measure map is isometric by definition. 
At this point one would expect to show the completeness of ([^x.c], D) in order to conclude 
the unique existence of a fixed flow of T by the Banach Fixed Point Theorem. This however 
becomes unnecessary, considering that we have an isometric isomorphism \x : [^x,c] Vx 
between flows and measures, which commutes with T by Theorem 3.5. 
By Theorem 2.3 and the remarks we have made on it, there exists a unique V G Vx '■ = v. 
By Theorem 3.5 there exists a unique [if)] G [^x,c] '■ nW\ — v and 

//[T0] = Tfi$\ = Tu = u = n$] T^ = 

since \i is bijective. Now let us suppose indirectly that [if)] is not a unique fixed point, 
meaning there is another if) ^ if) for which Tif> = if). Then 

/i[ip] = fj,[Tip] = Tfi[ip] => n[ip] =v = ii[if)] if)=if) 

which is a contradiction. 

Now let ipo G ^x,c be any flow, and consider the recursion ipi+i = Tif)L G *$?x,c and define 
vl '■= (J>[i>L\- Then vl+i = /^[T^l] = Tz/^ so by Theorem 2.3 we have that 31im(i/i) = z7 in 
(Vx,d). Since D([ ■],[•]) = d(fx[ •],//[• ]) and \i is injective, we have that 31im([^i]) = [if)]. 
We go on to characterizing the exact form of if). From Section 2.5 we recall that 

17 = lim y w a T*5 p 

\a\=L 

From the proof of Theorem 3.4 we see that if) c G ^x,c '■ Aif> c = T c Sq so 

if) = if, c * v = lim V] w a ip c * (S p o T" 1 ) = lim V] w a ip c (- - T a (p)) □ 

\a\=L \a\=L 

We remark that Aif> = a.e. and if Aip = a.e. then Aip^ = a.e. (L EN). 
Furthermore, taking if) := ip> CtP 

T^cji = ^2 W a T a^CP = ^2 W ^(' ~ T «G°)) ( m ° d 1c) 
\a\=L \a\=L 
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The first equality is clear, while the second equivalence requires some consideration. It can be 
shown by induction that z — T a {p) = ip a {T~ l (z) —p). So since if) CtP = ip c (- —p) we have 

V^pC^O)) = M T a\ z ) ~P) = M z ~ T a(p)) - (mod q c ) 

Summing over all addresses \a\ = L with factors w a we get the earlier congruence. 

Examining the intriguing quantities ipc{fa) further, using the notations of Definition 2.5, we 
see that 

w a ip c (ip a ) = w a i9 a + — - w a ln\ a (mod q c ) 

\a\=L \a\=L \a\=L 

Choosing the weights w k = \ s k (where s > 0, A| = 1 is the so-called pseudodimension 
of the IFS fractal (Ti, . . . ,T n ) which corresponds to its Hausdorff dimension if the open set 
condition holds) we have as the factor of ^ above the following quantity 

\a\=L \a\=L 

Considering the weights w k = X s k to be probabilities, this is the Gibbs Entropy Formula for 
a collection of classical particles. This can be interpreted as the entropy of the set of points 
F L = {T a {p) : a 6 A, \a\ = L} at each level L e N during the Chaos Game evolution towards 
the IFS fractal attractor F = (T 1; . . . ,T n ) = Cl{T a (p) : a e ^4}. Remarkably - plays the 
role of Boltzmann's constant, hinting at some potentially deeper interpretations. 

3.6 The Evolution towards Invariance 

As we have reasoned in Section 3.2, the intermittent evolution towards a fully developed 
turbulent flow field, corresponds to the iteration of the weighted / probabilistic transfer 
operator towards the invariant flow ip characterized in Theorem 3.6. 

One of the advantages of our model over classical Chaotic Advection formalism, is that the 
intermittent / iterative steps need not occur at equal time intervals, and in fact the spacing 
in time can be arbitrary. An interesting question is whether the fully developed state can 
be achieved in finite time. The answer is, it certainly can. If the time spacing of each 
iteration decreases say geometrically, then the total time to reach the L = oo level iteration 
(corresponding to ip) will be finite. This is clearly true for any convergent sequence of time 
spacing. 

Another question which may arise when considering the previous sections, is how the global 
superposition of n pushforwarded stream functions according to the transfer operator T re- 
lates to the local experimental picture of the splitting of eddies. Upon some contemplation, 
we may realize that the global and local viewpoints are the direct consequence of the as- 
sociativity of transfer T(T L ^ C)P ) = T L (T^ Ci p). Executing an iteration of T over the entire 
flow field T ip c ,p corresponds to the local intermittent splitting of the initial eddy ip CtP into 
n eddies Tip CtP and therefore the splitting of all of its level L iterates. 



23 



3.7 Equilibrium Points of the Invariant Flow 

3.7.1 Sink Singularities 

From the representation shown for the invariant flow 

if)[z) = lim w a ip c (z - T a (p)) 

L->oo ' * 
\a\=L 

where the limit is taken in the sense of ([^x,c]j D), we see that ip is an infinite shifted 
superposition of the eddy ipc of character c, preserved by if). Since ip c is centered at the 
origin, each weighted eddy w a i/j c (z — T a (p)) is centered at T a (p). 

Therefore the set of sink singularities of ip is precisely the countably infinite IFS fractal 
(Ti, . . . , T n ) = {T a (p) : a £ .4.} where p is the primary fixed point. Since as we have reasoned 
ip represents a fully developed turbulent flow field, this resolves Mandelbrot's Conjecture 2.1. 
Note that IFS fractals mostly have a non-integer Hausdorff dimension, but in certain special 
cases, such as Dragon Curves, the dimension can equal two. 

Since the flow concentrates on an IFS fractal, we must focus our efforts on the study of such 
sets, in order to uncover the characteristics of fractal potential flows. We also emphasize 
that by Theorems 3.1 and 3.2, having similarity contractions is necessary for preserving the 
physicality of flows under transfer. 

3.7.2 Saddle Points 

When visualizing T L if) C}P for large enough L £ N, one observes a thinning of the basins of 
attraction of each eddy center F L = {T a (p) : a £ A, \a\ = L} which raises the question 
whether the thinning continues on to a width of zero. The basins are partitioned by di- 
rectrices (manifolds) of the saddle points (hyperbolic equilibria) of the flow field. In the 
neighbourhood of each partitioning infinite separatrix (stable manifold) belonging to these 
saddles (the unstable one ending in sinks), the flow behaviour becomes chaotic - meaning a 
tracer particle starting down the streamline near one side of the separatrix, may end up at a 
distant sink singularity, relative to if it had started near the other side. When the flow field 
is considered as the phase portrait of a Hamiltonian system, then this signifies sensitivity to 
initial conditions near the particular separatrix. For the entire flow field to be considered 
chaotic, we must show that such partitioning separatrices "cover" the flow field, in some 
sense. We also intuit that this covering would be implied if the saddle points could be shown 
to be dense in the countably infinite IFS fractal (T 1; . . . , T n ). 

Conjecture 3.1 The infinite separatrices of the saddle points of ip are dense in the plane. 
Conjecture 3.2 The saddle points of ip are dense in (T 1; . . . , T„). 

We further conjecture that the second conjecture implies the first, which in turn guarantees 
the chaotic nature of the invariant flow field. 
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4 Visualization 



When considering the transfer operator over stream functions, one encounters a proliferation 
of branch cuts, as the transfer iteration progresses. This can make the visualization of 
streamlines - the level curves of the stream function - quite difficult. Note that stream 
functions in [^/x,c] are only required to be almost everywhere smooth, and the transfer 
iteration results in a countably infinite Hutchinson union of branch cuts in ijj, which pose no 
issue theoretically, despite the visual and algorithmic mess they create in practice. This is 
illustrated by Figure 2 with transfer parameters 



27T ; 2TT; 



Pi = 0, p 2 = 1, ipi = 0.65e s \ ip 2 = 0.65e 4 \ Wx >2 = 0.5 

Therefore in order to iteratively visualize the flow field, one seems to have no option but to 
resort to an iteration over the velocity field, ipo can be chosen arbitrarily, as apparent from 
Theorem 3.6, so we choose := J2k w kTk4> Ck , Pk as discussed at the end of Section 3.4, where 
Ck = Log (ft (choosing C = 2ir in Section 3.1) for an IFS of Section 2.4. The evolution in 
the velocity field progresses according to Theorem 3.2. So we have the following recursive 
iteration for the velocity field, over the flow space ([^/x,c]) D) with preserved character c = 
J2 k WfcLog V?fc and X D (T 1; . . . , T n ). 



vq{z) = v[ij)o] = ^w fe (Log if k )- 



Pk 



k=l 



Pk 



2 



vl+i(z) = v[T^ L ](z) = J2 Wk V~ la VL ^ T k 1 W) 

k=i mi 

Once we have arrived at a large enough iteration level L G N, we may execute a streamline 
solver algorithm, resulting in a preferable image as in Figure 3. The transfer iteration is 
visualized in Figure 4. Note the apparent convergence to an attractor flow field, as predicted 
by Section 3. 

In order to visualize the level curves of the corresponding potential function - the equipoten- 
tial lines - we first note the correspondence in character char(^) = — i char(^) by Definition 
3.1. If one does not wish to deal with the arising sources in the conjugate to the above ipo (as 
we did in Figure 5), we may choose simply instead ipo := ip CtP and thus ipo = ip-i c , P with the 
primary fixed point p and the above c, as discussed in Section 3.5. Then in order to arrive 
at the potential field at level L, we execute the same velocity recursion on v (z) := v[ip ] as 
above on Vq(z) = v[ipo]. The iterative formula remains the same because of Theorems 3.1 
and 3.2 implying the following, with v L := v[iPl), L G {0} U N. 

~P 



v (z) = v[tpo)(z) = -ic 



2 



VL+1 



ufe+i] = v[t^ l ] = v[Tip L } = 5^ Wfc rn2 T k v \^L\ = $^ Wfc r~i2 Tk * VL 



Depending on the sign of Im c, it may be more convenient from the viewpoint of a streamline 
solver, to take the negative of vq so that Re char(-u ) < 0. 
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Figure 5: The equipotential lines at iteration level L = 6, and a detail image. 
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5 Conclusive Remarks 



After attempting an extensive and analytical overview of Fully Developed Turbulence (FDT) 
in Section 2.1, we proceeded to introduce some theoretical concepts which we fused together 
in Section 3 in the form of Fractal Potential Flows. Our goal had been two-fold: (1) to formu- 
late the random intermittent evolution towards FDT as the iteration of a weighted transfer 
operator; (2) to show that this operator has a unique attracting fixed point, modelling the 
FDT flow field. 

As far as suggestions for future research, we mention the conjectures stated in Section 3.7.2, 
as well as the inspirational remarks on entropy in Section 3.5. Furthermore, the isometric 
isomorphism introduced in Section 3.4 between flows and probability measures, may be ex- 
ploited to translate former results on measures to the language of flows, while examining the 
physical implications. Considering Section 3.7.1, the study of Fractal Potential Flows (and 
thus FDT) reduces essentially to the geometry of IFS Fractals with similarity contractions, 
so their investigation appears to be of fundamental importance. 

We wish to express our appreciation to Dr. Sommerer and Dr. Tel for their clarifications of 
their work with Dr. Ott and Dr. Karolyi et al. and its relation to Aref's theory of Chaotic 
Advection. We express gratitude for the support of Prof. Edward R. Vrscay whose unbiased 
skepticism also motivated the author to chisel this attractor and its literature foundations 
to an arguably sufficient level of iteration, with his never sharp enough IFS maps. 

We dedicate our paper to the memory of Profs. Benoit B. Mandelbrot and Hassan Aref. 
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